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ABSTRACT 

Expressions  were  derived  for  the  laesn  and  mean  square 
value  of  the  output  of  a  phase  detector  to  vhich  is  applied  a 
sinusoidal  signal  corrupted  by  narrov-band  Gaussian  noise, 
together  with  an  oscillator  reference  of  the  saiae  frequency  as 
the  signal  b\xt  displaced  in  phase.  The  analysis  wets  based  on 
the  nodel  of  a  bailanced  phase  detector  coiQiosed  of  peak-detecting 
diodes  producing  a  difference  output  froa  the  two  halves.  Based 
on  the  above  output  paraueters,  a  aeasure  of  phase  detector 
performance  was  defined  in  terms  of  a  “noise -to -noise  ratio," 
the  square  root  of  the  output  variance  without  signal  divided  by 
the  same  quantity  with  the  signal  present.  The  case  \diere  signal 
and  reference  are  in  quadrature  was  treated  in  detail,  yielding 
curves  of  this  ratio  under  varying  input  and  reference  conditions. 
In  addition,  the  voltage  oxitput  "signal-to-noise  ratio"  was 
calculated  as  a  function  of  the  phase  angle  between  signal  and 
reference,  and  the  resulting  graphs  show  the  effects  of  reference 
level  and  input  slgnal-to-noise  power  ratio. 
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1 

A 

1.  IlfTRODUCnOH 

With  the  increasing  use  of  phase -con^parison  techni<;[ues  in 

the  design  of  guidance  and  connunlcatlcn  equipment,  it  is  lagportant 

to  know  how  devices  providing  such  operation  behave  under  vsurylng 

conditions  of  the  signal  and  corrupting  noise  envlronaent.  Idealljr, 

a  phase  ”coaq>«u:ator”  or  "detector”  conibines  a  received  signal  and 

a  locally  generated  reference  waveform  in  some  manner  to  produce 

an  output  voltage  idiich  is  proportional  to  the  phase  difference 

between  these  inputs.  Bbwever,  in  practice  this  desired  linear 

operation  does  not  continue  indefinitely  as  the  phase  difference 

increases  in  magnitude ,  instead  resulting  in  the  familiar  leveling 

1  2 

off  and  subsequent  decrease  in  detector  output  voltage.  ^  In 
cases  where  the  received  and  reference  waveforms  are  relatively 
unperturbed  by  noise  and  differ  little  in  their  time  behavior, 
the  non-linear  portion  of  the  0iase-d-*iector  characteristic  can 
be  Ignored  and  a  linear  analysis  of  its  performance  caurrled  out.’' 

But  with  the  advent  of  systems  worklxig  in  a  hif^  xiolse  environment, 
large  fluctuations  are  superiiQosed  on  the  received  signal  which 
cause  corresponding  phase  deviations  from  the  normally  steady 
reference  waveform,  resulting  in  excursions  beyond  linear  detector 
operation.^ 

It  was  for  the  purpose  of  considering  the  performance  of 
phase  detectors  under  such  high  level  noise  conditions  that  the 
present  study  was  undertaken.  In  azialysls  previously  carried  out, 
the  siiQ>llfying  assumption  was  made  thut  the  rexerence  amplitude 
was  much  larger  than  the  received  signal  magnitude,  resulting  in 
the  phase  detector  outpxxt  voltage  being  independent  of  this 
reference  a>iQ>lltude.^  While  this  is  a  valid  approximation  in 
practice  for  relatively  unperturbed  signals,  it  no  longer  holds 
when  the  added  noise  fluctxiations  are  sufficient  to  make  the 
over-all  received  asq>litude  coiqparable  to  the  reference  magnitude. 

The  present  investigation  proceeded  imder  gaxieral  signsQ.,  noise, 
and  reference  parameter  conditions. 
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2.  MODKL  OF  PHASE  DETECTOR  OPERATIOH 


noise  perfoxnance  of  the  idiase  detector  is  st^idied 
here  on  the  basis  of  a  nodel  of  its  operation  which  has  been 
used  in  previous  applications  to  autoaatic  frequency  and  phase- 
control  systeas.^  It  is  assuaed  to  be  a  balanced  phase  detector 
consisting  of  two  peak-detecting  diodes  and  associated  filter 
circuits  j  as  shown  below  in  Figure  1,  together  with  the  usual 


FIGURE  1  -  SCEEKATIC  DIAGRAM  OF  PHASE  DETECTOR 

transfozner  windings  to  Introduce  the  received  and  reference 
voltages.  *  The  esipresslon  for  the  output  voltage  e^  of  the 
phflise  detector  can  readily  be  obtained  froa  a  vector  diagram  of 
the  circuit  voltages,  once  the  latter  are  written  in  teras  of 
aaplitude  and  phase. ^ 
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For  the  present  noise  analysis,  the  received  voltage 
appearing  within  the  diode  circuits  is  taken  to  consist  of  a 
cosinusoidal  signal  of  angular  frequency  oi,  upon  which  is  super- 
ia5>o8ed  narrow-hand  Gaussian  noise  centered  at  aj.  By  the  well- 
known  Rice-Bennett  representation,  e.  can  be  expressed  in  the 
desired  polar  form  by: 

e^  =  p  cos  9^  ,  (1) 

where  p  Is  the  azoplitude  of  signal  plus  noise: 

p  =  +  x(t)^^  +  y^(t)  ,  (2) 

A  =  the  amplitude  of  the  received  signed., 

3c(t);  y(t)  =  the  slowly  time-varying  in-phase  and 

quadrature  amplitudes,  respectively,  of  the 
noise  centered  at  cu,  idaich  are  independent 
Gaussian  random  variables  with  mean  value 
zero  and  standard  deviation  0,  and 

0^  =  the  tot«LL  noise  power  (for  unit  resistance); 

9^  =  cut  +  9^  ,  (3) 

9^  being  the  phase  of  signal  plus  noise: 


The  reference  voltage  e  within  the  diode  circui.ts  is  a 

r 

sine  wave,  assumed  to  have  the  same  angular  frequency  as  the  signal 
portion  of  e^  and  an  arbitrary  phase: 

=  B  sin  92  ,  (5) 

where  B  =  the  aa5>litude  of  the  reference  voltage. 
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9^  being  the  arbitrary  phase  angle.  Thus,  vhen  9^  =0,  is 
in  quadrature  with  the  signal  portion  of  e^,  a  condition  typical 
of  sany  phase  detector  applications. 

^ntie  filter  circuits  associated  with  the  diodes  (See  Figure  l) 

are  sade  with  sufficient  bandwidth  to  pass  the  narrow-band  noise 

coB^nentSj  and  at  the  saae  tlae  reproduce  the  crests  of  the 

half-wave  rectified  voltage.  !I5ius,  the  peak-detected  contributions 

of  both  halves  of  the  balanced  circuit  are  seen  froa  the  vector 

diagram  of  the  circuit  voltages  in  Figure  2,  where  their  positions 

6 

take  into  account  the  quadrature  relation.  To  one  of  the  diode 


FIGURE  2  -  VECTOR  DIAGRAM  OF  CIRCUIT  VOLTAGES 

circuits  is  applied  the  vector  sum  of  e^  and  half  of  the  reference 
voltage  ^JjZ  during  its  conduction  interval,  while  the  other  receives 
the  vector  difference.  Upon  peak  detection,  each  half  produces  a 
rectified  filtered  output  voltage  equal  to  the  magnitude  of  the  total 
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applied  voltage  vector.  Noting  the  magnitudes  of  the  individual  vectors 
on  Figure  2  (obtained  from  Equations  (l)  and  (^)),  and  defining  the 
phase  difference  (p  between  e^  and  e^  in  addition  to  that  of  quadrature: 

<P  -  <Pl  -  <P2  > 


then  the  output  voltages  e^  and  of  the  two  halves  are  obtained 
from  the  triangle  law: 


2 

p  +  +  Bp  sin  cp 

2 

p  +  - Bp  sin  Cp 


1/2 

1/2 


(8) 

(9) 


where  the  positive  square  root  is  always  understood^  and  from 
Equations  (3)#  (6),  and  (?): 


9  »  9^  +  9q  •  (^) 

Since  the  ladividual  peah  voltages  appear  across  the 
output  terminals  in  opposite  directions,  the  phase  detector 
Otttpttt  Tolttse  consists  of  their  aifference: 

=0  '  'l  ■  ®2  • 

Equations  (8)  through  (ll)  give  the  analytical  representation  of 
phase  detector  behavior-  to  be  used  in  subsequent  calculations . 
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;1.  STAIISTICAL  PR0P2SIIE3  OP  T3E  CXJTFJT 


Sie  expxsssiona  derived  for  tlie  phaae  detector  ahoif  that 
the  o\xt|mt  Toltflge  depends  on  the  statistical  nature  of  the  noise 
si^riaqposed  tqpon  the  xeceiTcd  signal.  As  a  result^  the  performance 
of  the  phase  detector  can  only  he  described  in  texms  of  probability 
distributions  and  their  associated  aoaents.  Tuo  of  the  most 
useful  factors  related  to  these  are  the  "asan"  and  "laean  square" 
values  cf  the  output,  yhich  are  treated  in  this  analysis  in  that 
order. 


Noting  froa  Equations  (8)  through  (U)  that  the  output 

voltage  e^  is  e:Q>ressed  in  tens  of  the  amplitude  p  and  phase 

of  the  signal  plus  noise,  it  is  necessary  to  knov  the  Joint 

probability  distribution  of  these  randcmi  variables  in  order  to 

fonaUate  the  desired  output  properties.  Such  a  density  function 

is  obtained  from  the  corresponding  description  of  the  in-phase  and 

q,uadrature  amplitude  variables  of  the  corrupting  narrov-band  noise, 

X  and  y,  respectively.  It  vas  assumed  that  they  vere  indepe!t3S!lt 

Oaussian  randon  variables  vith  mean  value  zero  and  standard 

deviation  (5,  uhere  is  the  sum  of  the  equal  individual  noise 

posers  (T/Z  of  the  two  components,  called  'ihe  "total  noise  pover." 

0 

!Ehus,  the  Joint  distribution  V2(x,y)  takes  the  simplified  form: 


W2(3£,y) 


(12) 


Sy  means  of  Equations  (2)  and  (4)  treated  as  transformation 
relations,  the  desired  Joint  distribution  of  amplitude  and 
phase  of  the  signal  plus  noise  is  obtained^: 


2  2 

(p  +  A  -  2Ap  cos  9^) 


(13) 
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The  application  of  Xq',i&tlon  (13)  and  the  expression  for  e 
ter»8  of  p  and  cp^  to  the  veil-known  statistical  averaging 
foraulas  yields  the  mean  value  and  mean  square  value 
of  the  phase  detector  output: 


o 


-  2« 


i 


o  Jo 


W2(p,9i)dPidp 


(11^) 


r 


JO  Jo 


2jt 

Je^(p,(p^)]^  W2(p,<p^)d9^dp 


(15) 


^ere  the  integrations  are  carried  out  over  the  appropriate  ranges 
of  the  aaq>litude  and  phase  random  variables. 

While  the  desired  calculations  can  be  carried  out  directly 
with  Xquatioxxs  (l4)  and  (15),  it  is  more  convenient  to  convert 
these  eapressions  back  into  the  original  in-phase  and  quadrature 
noise  variables.  As  a  result,  the  following  equivalent  pair  of 
formulas  is  obtained: 

CO 

eQ(x,y)  W2(x,y)dxdy  ,  (16) 


’  00 


r 


00 


W2(x,y)dxdy 


''-00  ■'-00 


(17) 


\diere  x  and  y  have  infinite  ranges,  and  the  phase  detector  output 
voltage  is  ex])ressed  in  terms  of  them  with  the  Gdld  of  Squatlons 
(2),  (4),  and  (8)  throiigh  (u): 


j 
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L' 


r 


B  \ ^  B 

X  +  A  +  ^  sin  ®  +  y  +  o  CO® 

cot  «  o; 


1/2 


B  i  B 

X  +  A  -  sin  9J  +  y  -  ^  cos  (pi 


.21 


1/2 


J 


(19) 


Together  with  Xquation  (12),  Equations  (16)  throu^  (18)  form  the 
baais  for  the  calculation  of  and  ^  in  this  study. 

3.1  THE  OUTFIT  HBAB  VALUE 

To  aid  in  the  application  of  Equation  (16),  it  was  found 
advantageous  to  express  the  phase  detector  output  voltage  e^(x^y) 
in  integral  fora.  The  desired  expression  resiU-ts  fro»  a  table  of 


Laplace  transforas: 


.9 


e  (x,y)  * 

j 


r*  \ 

sixih  <Bt  [(x-fA)8in  q»^+ 


y  cos  ^ 


(19) 


vhich  is  raUd  over  the  infinite  ranges  of  integration  for  x  and 
y.  Based  on  the  8d}ove  representation ^  the  detedled  evaluation 
of  ^e^^  is  carried  out  in  Appendix  A,  vhich  yields  the  foUovlng 
expression  for  the  output  aean  value: 

-  og)  (S^)+  (a-e)l^  (5^j 

Tdiere  I^Cx)  and  I^(x)  are  aodified  Bessel  functions  of  the  first 
kindj  and  a,^  are  the  noraalized  paraaeters : 


(20) 


1 

26^ 


.2  B 

+  ip  ; 


(21) 
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ExajBination  of  Ekiuation  (20)  with  the  aid  of  Equati  a  (2l) 
shows  that  ^e^^  has  the  expected  behavior  for  extreiae  values  of 
the  above  paraswters,  based  on  the  model  xised.  Thus,  when 
either  the  signal  or  reference  aarplitude  vanishes,  or  these  two 
sinusoidal  caB5>onents  are  in  quadrature  (9^  =  O),  the  output 
mean  value  is  zero.  On  the  other  hand,  in  the  absence  of  noise 
(6  =  0),  the  phase  detector  yields  a  mean  value: 


which  for  the  case  ^en  signal  and  reference  are  in  phase 
(9q  =  k/2),  reduces  to  just  the  referesnce  amplitude  B. 

3.2  THE  OUTPUT  MEAN  SQUARE  VAIUE 


As  in  the  case  of  the  mean  value,  the  square  of  the  phase 
detector  output  voltage  was  expressed  in  integral  form  to 
facilitate  the  evaluation  of  Equation  (1?).  This  is 

obtained  frcsa  the  following  equivalent  expression  for  e^^(x,y), 
derived  from  Equations  (8)  through  (U): 


eQ(x,y) 


2Bp  sin  (9^  +  9^) 


(33) 


From  the  table  of  Laplace  transforms,  the  resulting  Integral  is: 


.9 


eQ(x,y)  =  2B  [(x  +  A)  sin  9^+7  cos  9^j  ♦ 

00 

j  I 

~  /st  [(x+A)  sin  9^  +  y  cos  9^]  \  • 

^  o  '  ’ 

•  exp  t  [(x+A)^  +  ^  ^ 

which  is  again  valid  over  the  infinite  ranges  of  integration  for 
X  and  y.  As  derived  in  Appendix  B,  the  resulting  expression  for 


(2h) 
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is  ®  coHplicated  pover  series  in  jsee  Equation  (2l)] ,  >diose 
coefficients  contain  integrals  involving  confl^^ent  hypergeoartric 
funct.^on8;  for  idiich  no  exact  closed  form  evaluation  can  be 
found.  Consideration  of  this  general  foxmila  viU  not  be 
undertaken  here,  since  of  prime  Interest  are  calculations  of 
phase  detector  performance  vhich  are  treated  in  the  following 
sections . 
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4.  CALCULATICJN  OF  Oin?PUT  HOISK -TO -NOISE  RATIO 

While  the  mean  and  mean  sqviare  values  previously  detemined 

give  a  partial  statistical  description  of  the  phase  detector  output, 

these  factors  in  themselves  do  not  provide  a  criterion  for  the 

performance  of  sxich  a  device  in  the  presence  of  noise.  In  general., 

a  detection  system  has  two  states  of  interest,  one  vhen  a  signal 

is  not  present  and  one  vhen  a  signal  is  present*  From  this  it 

is  natural  to  seek  some  parameter  with  vhich  to  coiqpare  the  effect 

of  the  ever-present  noise  on  both  states.  Since  in  each  case 

there  are  output  fluctuations  about  some  mean  value  (vhich  may  be 

zero),  a  sxiitable  choice  vould  be  the  "variance  of  the  output," 

2 

denoted  by  .  In  terms  of  the  moments  previously  discxissed,  this 
quantity  is  defined  by: 


For  the  case  ^ere  the  signal  is  present,  the  expressions 

derived  in  Appendices  A  and  B  ere  sxibstituted  directly  into 

2 

Equation  {Z'y) ,  and  the  resxilting  variance  is  denoted  by  0  . 

s 

When  the  signal  is  not  present,  an  evaluation  of  these  formulas 

at  A  =  0  must  first  be  made,  with  the  variance  thus  obtained 

2 

being  denoted  by  (5  ,  since  only  noise  is  received  here.  From 

these  tvo  definitions,  the  criterion  of  performance  for  the  phase 

detector  is  taken  to  be  the  noise-to-slgnal  plus  noise  output 

variance  ratio,  to  the  one -half  pover.  This  quantity, 

10  11  “  ° 

is  called  more  briefly  the  "noise -to -noise  ratio."  '  Since 
the  output  nKian  value  vas  seen  to  be  zero  for  A  =  0,  this  ratio 
can  be  vritten  in  terms  of  the  moments  from  Equation  (25)  in  the 
following  form: 
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vhere  the  evaliiatlon  at  zero  signal  asqplltude  Is  indicated  above. 
Equation  (26)  is  now  applied  to  an  operating  case  of  interest 
for  the  phase  detector^  nanely  ^en  the  signal  and  reference 
are  In  quadratxure,  i.e.  9^  =  0.  This  is  additionally  desirable 
because  the  output  signal-to-nolse  ratio  is  not  meaningful  vhere 
no  measure  of  signal  exists. 

In  this  case  of  null  operation  in  the  absence  of  noise, 
the  output  mean  value  (e^  =  0,  even  with  signal  present.  From 

Equation  (21),  p  =  0,  so  that  from  Equation  (Blfl)  of  Appendix  B 
the  expression  for  reduces  to  the  single  integral: 


\ihere  stands  for  the  confluent  hypergeometric  Itmction,  a  is 
defined  in  Equation  (2l),  and  r  is  the  diaensionless  ratio: 


(28) 


Setting  A  =  0  in  Equation  (2?)  and  substituting  it  and  the  original 
form  into  Equation  (26)  with  =  0  yields  an  exact  expression 
for  the  output  noise -to-noise  ratio. 


To  relate  the  calculations  s»de  with  Equation  (26)  for 
this  quadrature  case  to  laboratory  laeasurements  on  actual  devices, 
the  above  parameters  are  eacpressed  instead  in  terms  of  the  total 
received  power  P^,  the  reference  power  P^,  the  input  signal~to -noise 
power  ratio  ”a,”  and  the  reference -to-total  received  power  ratio 
”b".  These  are  given  by  the  relations  (for  unit  resistance): 
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As  a  resill,t,  t^e  paraaieters  a  aod  r  appeeuring  in  Equation  (2?) 
take  the  fom: 


s  ^  (1  +  s)  } 

(31) 

b(l  +  a) 

(32) 

The  evaluation  of  the  output  noise-to-noise  ratio  using 

the  exact  integral  foraula  of  Equation  (2?)  led -to  introlved  tine- 

consuaing  nuaerical  procedures >  so  instead  various  approximation 

methods  were  considered^  the  details  of  ^ich  appear  in  Appendix  C> 

A  si3i5>lified  functional  fom  which  fitted  the  confluent  hypergeonetric 

function  over  a  rax2ge  of  the  parameter  r  led  to  an  fiqppro>JLaate 

12  i 

result  in  terms  of  the  modified  e:q>onential  integral.  But  since 
the  latter's  numerical  evaluation  hecsaae  quite  Involved^  further 
sii5>lifications  were  sought.  As  a  result,  it  vss  possible  to 
o)>tain  reasonable  closed-fozm  e:q>re88ion8  for  two  ssysQxfcotic  cases. 
These  are: 

Reference  Amplitude  Halt  (b — )0): 


(33) 


Leuge  Reference  Amplitude; 


n 


1  _  1  2a  1$  -  9(3-wt)  +  ^{6■*6ei^^er) 


2r 


2r 


1/2 


(3i^) 


where  r  is  expressed  in  terns  of  "a”  and  "b"  by  Equation  (32). 
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The  preceding  fonatLlas  are  presented  graphically  in 

Plgure  3  ia  the  fora  of  curves  of  6  /(S  for  the  quadratvire  case 

(<Pq  *  0)  as  a  function  of  "a”  for  tj.e  lower  limit  and  large 

levels  of  "b”.  Intermediate  reference  values  would  result  in  plots 

falling  in  between  the  above  extremes,  and  these  would  be  appropriate 

for  subsequent  comparison  with  experimental  data.  But  the  curves 

calculated  here  do  reveal  the  principal  effect  of  the  signal  on 

the  phase  detector  in  decreasing  the  output  noise  over  that  in 

its  absence,  as  seen  in  Figure  3  by  the  6/6.  V£a.ue8  staying  above 

n  e 

unity..  This  is  (jposite  to  the'hehavior  of  aarplitude  detectors, 
so  judging  from  the  output  noise-to -noise  ratio  it  can  be  concluded 
that  the  phase  detector  does  not  operate  in  the  manner  normally 
ascribed  to  detection  processes. 

The  spread  between  the  curves  of  Figure  3  illustrates  t.. 

relative  merits  of  the  two  extreme  modes  of  operation.  For  large 

reference  levels,  the  6j6^  ratio  becomes  insensitive  to  input 

ns 

slsnal-to-noise  ratio  changes,  but  never  rises  much  above  unity. 

On  the  other  hand,  the  use  of  a  low-reference  level  results  in 
significant  oittput  noise  decreases  with  the  signal  present,  wliich 
could  serve  as  the  basis  for  scrae  t;ype  of  detection  scheme,  but 
here  the  6^/6^  ratio  is  affected  more  by  varying  input  conditions. 
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5.  CALCULATION  OF  OUTPOT  SIGKAL-TO-NOISE  RATIO 


For  chai*acterlziDg  the  response  of  a  physical  device  to 

a  conibination  of  signal  and.  noise  waveforms^  the  laost  coiaaionly 

used,  criterion  of  perforaance  is  the  "output  signal-to -noise  power 

ratio,"  denoted,  by  "a^"*  For  the  type  of  input  to  the  phase 

detector  and  with  the^ output  specified  in  terias  of  voltage,  the 

lU  15 

square  root  of  the  above  quantity  is  thus  defined  as;  ^ 


average  output  for  signal!  raverage  output  for 
plus  nol^  input  J  ~  L  noise  input  only 
rstaniard  Aviation  of  the  output  acout!. 

L  the  nean  for  signal  plui>  noise  input  J 


Since  froa  Equation  (20)  it  is  seen  that  for  no  input  signal  present, 
i.e.  A  =  0,  =»  0,  then  In  terns  of  the  nonents  previously- 

discussed,  this  ratio  is  gi-ven  by: 


As  indicated  pre-vriously,  -the  exact  evaluation  of  do*8 

not  have  a  fora  .suitable  for  application  to  Equation  (36),  but  more 
tractable  approximate  e:q>ressions  were  ob-tained  for  limiting  cases 
of  -the  reference  paraae-ter  "b."  Specifically,  the  derivations  of 
Appendix  -D  yielded  for  the  -voltage  output  slgnal-to-noise  ratio 
in  -the  general  case  of  an  arbitrary  phase  angle  separating  -the 
signal-reference  quadrature  condition; 


(37) 
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Small  Beference  Aaplltude  Limit  (b — >0): 


^ere  "a"  is  the  input  signal-to-noise  power  ratio  defined  in 
Equation  (30),  and  and  are  modified  Bessel  functions  of  the 
first  kind. 

Graphs  of  the  calculations  with  Equations  (37)  and  (38)  are 
given  in  Figure  4  in  the  form  of  curves  of-^^  as  a  function  of 
the  sine  of  the  phase  angle  for  two  values  of  ’’a"  and  the  limiting 
reference  levels.  In  spite  of  the  radically  different  functional 
forms  exhibited  hy  the  preceding  two  expressions,  their  resulting  plots 
are  seen  to  approach  each  other  very  closely  over  the  Initial  angle 
range  for  each  va3v.e  of  input  signal-to-noise  power  ratio.  Thus, 
such  a  relatively  narrow  spread  shown  in  Figure  4  would  appear  to 
include  the  curves  at  intermediate  reference  levels,  which  should  he 
the  subject  of  coaputational  verification  in  the  future. 

For  the  first  6°  off  the  signal-reference  quadrature  condition, 
the  voltage  output  signal-to-noise  ratio  varies  linearly  with  sin  <5)^, 
with  the  constant  of  proportionality  being  *^Ea  or  subtly  below 
over  a  wide  range  of  reference  levels.  With  larger  phase  offsets, 
the  USmitlng  curves  will  diverge  more  markedly  due  to  the  subsequent 
departure  from  linearity  of  Equation  (38)*  But  with  the  phase  detector 
t^yplcaily  operating  in  qiiadrature,  a  small  phase  shift  occurring 
between  signal  and  reference  will  produce  an  increasing  output  signal-to- 
noise  ratio  wtxich  provides  the  means  for  cor. -cting  such  waveform 
displacements . 


6.  CONCLUSION 
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The  analysis  of  phase  detector  performance  for  a  signal 
plus  noise  input  was  carried  out  for  a  halaaced  circuit  with  two 
peak-detecting  diodes.  The  model  of  narrow-hand  Gaxissian  noise 
used  was  considered  appropriate  for  most  of  the  applications 
encountered.  The  criteria  of  "noise -to-noise  ratio"  aid 
"signal -to -noise  ratio"  at  the  output  were  selected  to  evaluate 
the  utility  of  this  detection  process.  Although  computational 
difficulties  precluded  obtaining  a  con^jlete  range  of  analytical 
resvxlts,  the  approximate  formulas  which  were  four.d  served  to 
illustrate  the  salient  features  of  phase  detector  behavior. 

For  the  signal  and  reference  in  quadrature,  a  relatively 
small  decrease  in  output  noise  results  when  the  signal  is  present 
for  a  wide  range  of  input  signal-to -noise  power  ratios  at  large 
reference  levels,  but  a  larger  spread  of  output  noise-to-noise 
ratios  is  exhibited  at  low  reference  levels  under  varying  input 
conditions.  When  the  signal  and  reference  depart  somewhat  from 
the  quadrature  condition,  the  voltage  output  signal-to-aoise  ratio 
is  roughly  proportional  to  both  the  sine  of  this  phase  displacement 
and  the  corresponding  input  ratio  (for  valuec-  between  lUiity  and  ten) 
over  a  wide  range  of  reference  levels. 

It  is  the  aim  of  this  scaacu_^fc  restricted  analysis  to 
encourage  further  efforts  in  investigating  more  general  phase 
detector  performance  both  from  a  thsoreticeQ.  and  experiaential 
viewpoint.  One  ouch  study  has  considered  the  same  t^/pe  of  circuit 
with  noise  present  in  both  the  signal  end  reference  channels,'*'  but 
the  results  are  too  coaiplicated  to  yield  practical  ii::fox*3jatlon  and 
in  eiddition  stress  the  square -law  behavior  of  the  diodes  which  is 
practice  is  a  less  realistic  nsodel.  Since  equipment  applieatiorvS 
are  tending  increasingly  towai'd  the  use  of  more  sophisticated  phare 
detectors  to  improve  performance  from  the  uoise  and  other  standpoints, 
it  becomes  essential  to  extend  the  presently  'snown  analytical  teihni,qxK8 
in  order  to  understand  the  behavior  of  such  devices . 
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APPENDIX  A 

-  EVALUAnOK  OP  THE  OUTPUT  MBAH'VAUJE 


The  integral  representation  for  the  phase  detector  output 
voltage  e^Cx.y),  Equation  (l9)>  and  Equation  (12)  for  the  Joint 
distribution  V^Cx^y)  are  sulsstituted  into  Equation  (l5)  and  the 
orders  of  integration  reversed.  By  the  addition  formula  for  the 
hyperbolic  sine: 


sinh(u  +  v)  =  slnh  u  cosh  v  +  cosh  u  sinh  v  ,  (Al) 


the  expression  for  /e  \  contains  only  products  of  integrals  in 

\  ®/ 

X  and  y: 


The  second  double  intsgrcd.  in  the  bracket  vanishes 
because  of  the  integration  with  respect  tc  y  of  an  odd  function 
over  synraetrical  limits.  If,  in  the  remaining  part,  tha  variable  x 
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is  translated  by  an  aaount  -  A  (the  infinite  limits  remaining  unchanged) 
and  the  square  is  completed  in  the  resulting  e^qionent: 


exp 


•  K* 


ll  j  a2 

1  L  1  ' 

.-4 

1  ) 

(A3) 


exp{- 


t  + 


26 


X  - 


t  + 


26^  I  2C^ 


-1- 


then  there  results  for  /e 


\  - 

1  ,  I 

dt  1 

b2  i 

-  iptj  exp 

1  A^ 

1  1  1+  4.  1  1 

V 

0  ^ 

[  26^ 

[  zA  2^\ 

1 

x(t)y(t)  .  (A4) 


where 


00 


X(t)  =1  exp< 


Y(t) 


-y  joi  '' 


t  + 


2(? 


X  - 


A  1  ' 

?l'  • 


sinh(Btx  sin  dx 


t  + 


26 


y  )  cosh  (Bty  cos  9^)  dy 


{A6) 


(A5) 


The  above  integrals  are  evaluated  by  Reference  17^  P*  l64, 

nos.  (3)  and  (4),  yielding: 

-ll 


X(t)  =  rr 


1/2 


^2. 2  .  2„ 

B  t  sin  9 


26^ 


.  ^  V  jABt 

•  sinh  ( — 5- 

26 


[' 


exp  \ - 1;- 

1-1 


t  + 


1 1 


?)  (■ 


y(t) 


n 


1/2 


■1/2 


t  + 


I  -02*2  2„ 

I  B  t  cos  9^ 


26 


exp 


-1 


t  +  -- 


26 


\ 

!  ■ 


(A7) 


(AS) 
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Multiplying  Equation*  (A?)  and  (AS)  together,  coabining  and 
sioplifylng  the  exponents,  the  following  expression  is  obtained: 


X(t)Y(t)  =  * 


t 


Zf? 


-1 


«*P  {--IT 

-1 


t  + 


2(5^ 


-l| 


(A9) 


*  sinh  t  +  ^ 


\Z(^  \  Zf? 


sin  9 


Noting  from  Equations  (Alf)  and  (A^)  that  in  the  eaqponents 

-1 


1--^ 

1  1"^ 

t  +  =  t 

t 

2(f 

2(5^, 

the  8\ibstitution  of  X(t)Y(t)  into  canceling  terms  yields 

the  aore  coopact  fora: 


<*o> 


dt 


-1 


e^q?  (-  Qtt  t  + 


26^) 


-i! 


^slnh  (pt 


t  + 


2(5^ 


,(Ali)) 


Mhere: 


.2 

A  +ip 


(All) 


The  evaluation  of  Equation  (AlO)  proceeds  from  Equation  (A9) 
by  making  the  substitution: 


(A12) 


t  +  -^1  -  26^  1 

1  -  t 

t  +  -ijj 

-1-. 

21^1  1 

. 

2S®1 

which  resiOts  in  as  the  difference  of  two  integrals; 
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ko)  -  ^[j 


sinh  /pt  t  + 


(A13) 


dt  .  ^  1 


-  at  t  +  — r  )  sinh  (pt  t  + 


ii-H 

2(5^  I 


In  the  fir&t  integral  an  integration  by  parts  is  perfonaed,  lfc*-ting: 


u  =  exp  { -at 


Taking  the  indicated  differential  and  integral  yields: 


du  =  — 9  exp  < -  at  t  + 
26  !  \ 


.  t  + 


oosh/pt  t  +  -~  I  - 

2(r  I  L  I  '  26*^^  f 


aBlnh  /pt  t  +  — X 


V  =  - 


' '  ■«  "If  * 


Substituting  into  the  parts  formula,  the  evaluated  part  vanishes  at 

-1/2 

t  =  «  due  to  the  t  '  factor  and  at  t  =  0  due  to  the  h^erbolic 
sine  factor  being  of  order  t.  Puttying  the  remaining  integral  part 
into  Equation  (AI3)  and  arranging,  the  result  is  expressed 
in  three  terras: 

<'fo)  -  f  ^ 


fz[ 


l_l* 

1 

( 


sinh  {pt  t  +  -~ 
i  '  26 


(Al4) 


t  + 

262| 


t  +  — exp  /  -  at  t  +  -~  N  sinh  <Pt  t  +  -^ 


26^'  I 
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Next,  the  foUovlxig  change  of  A'ailahle  is  nade; 

,  \-l 


u 


t  + 


=  dt 


t  + 


-2 


a  .  1  .  Viad-u)^^  .  (t  +  -i-1  .  3-  ..  ; 

2e^(i-u)  ’  ’  \  28^1  2a^(iru) 


at  t  *  0,  u  B  0;  at  t  B  «•,  u  *  1. 

Suhstltutlng  into  Equation  (Al4)  results  In  sloq>ler  Integral  forms: 

r  1 


VI 


2P 
L  J 

1 


e“<*u  ^ 


u 


-  aa 


u 


w 


sinh  pu  du 


uV2(i.u)V2 


€'’®‘  Sinh  Pu  du 


(AI5) 


Two  basic  Integrals  are  to  be  evjsluated  in  Equation  (AI5), 
one  being  defined  by  R(p): 

1 


R(p)  » 


r-PU 


du 


(AI6) 


This  is  foimd  directly  in  Reference  table  P*  138>  no.(l4), 
for  t  »  u,  b  »  1/2,  V  =  0.  Noting  that  n(3./2)  =  y~«,  there  resxilts 
the  evalU/ition: 


R(p)  =  , 


(AI7) 
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^er«  I^(p/2)  is  the  modified  Bessel  function  of  the  first  kind  and 
order  zero.  Next,  the  second  basic  integral  to  be  evaluated  is 
S(p),  defined  by: 


S(p)  =  I 


du 


(AlS) 


u 


First,  differentiating  under  the  integreLl  sign  with  respect  to  p  yields: 


^  .  u)V2  C-: 


■ptt 


dp 


du 


(A19) 


Applying  Reference  9,  table  It. 3,  p.  138,  No.  (14),  for  t  =  u, 
b  *  1/2,  V  =  1^  noting  that  T-.S/s)  =  there  results  the  evaluation: 


^isl  ^  .  ji  €-p/2  I  {£ 

dp  2p  -^ll  2 


(A20) 


vhe.re  I^(p/2)  is  the  modified  Bessel  function  of  first  kind  and  order. 
Equation  (A20)  can  be  integrrted  directly  by  using  Reference  18,  p.  57, 
3*105^  for  V  =  -1,‘u  =  p/2,  noting  that  I_^(p/2)  =  I^(p/2): 


kll'i  ^1(1 


+  c 


(A21) 


The  arbitrary  constant  is  eva^-uated  at  p  =  0,  yielding: 


C  =  S(0)  -  I 


(A22) 


tdiere: 


^  ^  / 

I  ..^l/2 


s(o)  =  I 


J 


w 


du 


(A23) 
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From  Reference  19^  p.  196,  no.  855‘1,  x  =  u,  m  =  1/2,  n  »  3/2, 
substituting  in  for  the  required  gais&a  functions: 

ril!  =  >  ril)  =1  .  r(3)  =  X  , 

S(o)  is  found  to  have  the  value: 


S(0) 


K 

2 


(A24) 


Substituting  Equations  (A22)  and  {AZk)  into  (A2l}  thus  gives  the 
coDQilete  evaluation  of  S(p): 


S(p) 


(A25) 


By  expressing  the  hyperbolic  sine  and  cosine  in  terms  of 
eixponentials  and  combining  ttiem  with  the  other  exponent,  Equation  (AI5) 
can  be  written  In  terns  of  the  basic  integrals  Just  evaluated  in  the 
following  rearranged  form; 


(a  +  ^)s(a  +  p)  -  (a  -  p)  s(a  -  3) 


+  I  R(a  +  3)  -  I  R(a  -  3) 


(A26) 


Substituting  Equations  (AIT)  and  (A25)  into  Equation  (A26)  with  the 
proper  arguments  and  collecting  terms,  the  evaluation  of  the  mean 
value  ^e^  of  the  phase  detector  output  is  coaplete: 

-  exp  +  {o^)iJS^) 

This  resixlt  appears  as  Eqxiation  (20)  of  tfee  xaain  text. 


APPENDIX  B 
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EVALUATION  OF  THE  OUTPUT  MEAN  SQUARE  VALUE 


The  integral  representation  for  the  square  of  the  phase 
detector  output  voltage  (x>y),  Equation  and  Equation  (l2) 

for  the  Joint  distribution  V^{x,y)  are  substituted  into  Equation  (IT) 
and  the  order  of  integration  reversed.  In  this  case,  however, 
products  of  integrals  in  x  and  y  are  not  immediately  obtained 
because  the  modified  Bessel  function  does  not  have  the  simple 
addition  formula  previously  seen  for  the  hyperbolic  sine  in  Appendix 
A.  Nevertheless,  a  separation  of  integrands  can  be  effected  if 
an  integral  representation  for  I^(z)  is  applied  from  Reference  20, 
p.  202,  no.  179>  for  n  =  1  and  z  =  the  given  argument: 


r* 


(Bl) 


a  .  i  !  d6  cos  Q  exp  Bt  cos  ©[(x+A)  oin  cp^+  y  cos  9  . 

o 

Putting  Equation  (Bl)  into  and  again  reversing  the  order  of 

integration,  the  following  product  integral  form  is  obtained: 


B 


r* 


dt 


B" 

F 


d©  cos  9  • 


(B2) 


'o 

CO 


sin  (p 


/  r 

t 

(x+A)  exp  /  - 


■  f 


t(x4A)^+  ^ 


2^ 


-  Bt  cos  0(x+A)sin  <p^\  dx 


00 


exp  (- 


X 


t  + 


2(5^ 


y  -  Bt  cos  6(y)  cos  <P  >  dy 


+  cos  <p 


exp  (- 


UOO 


t(x+A)^  +  ^ 

zcr 


0( 


-  Bt  cos  9(x+A)  sin  <p^>  dx 


Bt  cos  0(y)  cos  dy 
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If  in  tlxe  iategratlcns  with  respect  to  x.  this  variable 
ia  traaslated  by  an  aaount  -  A  (the  infinite  liMlts  1‘ettaining  tmchanged) 
and  teras  in  pcwera  of  x  are  collected  in  the  integrand,  then  there 
results  for 

L  2\  B 


7 


dt 

b"  1  1 

-f  exp 

0 

-r*j  J 

d9  cos  $ 


(B3) 


exp 


/  *2  1  / 


where: 


X2(e,t) 


Yj^(«,t) 


r 

U 

r 


1-1 

|t  4--^ 

1' 

1  26^/ 

+  I  Ajr  -  Bt  cos  B  sin  x )  dx 


X  exp  (- 


t  f 


26^ 


x2. 


A  ' 

--  -  Bt  coB^fiin  9  |x>  dx 


exp  <- 


J 


t  +  y^  -  (3t  cos  &  COB  9„)y)  dy  , 

2(r  ° 


(B4) 


(B5) 


(B6) 


CO 


i-l 


y  exp  < -  t  + 


^  ’  y^  -  (Bt  cos  ©cos  9  )yldy  . 


26"/ 


’0'''|' 
) 


(B7) 


The  above  Integra]^  are  cveiluated  by  Reference  17,  314,  p.  65,  nos. (6a)  and 
(6c),  suppressing  the  exp  (-c)  factor  on  both  sides; 


^  00 
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both  integrcds  beir<g  valid  for  a>  0;i  fhich  is  satis  if  ed  in  this 
derivation  for 


a  =  t  +  . 

26 


(BIO) 


Making  the  appropriate  autstitutiOTus  fox*  the  constant  b  froo 
Equations  (B^+)  through  (BT),,  there  resu3.t  the  evaluations: 


2/0  t  1  III  T  ^ 

'  t  -i-  exp  t  +  (  "o  ~  Bt  COB  S  sin  qs  )  , 

26^1  ri  26‘^  \(^ 

i  I 


(BU) 


“  I  Sir 


^  |Bt  cos  6  sin  f  -  ~ 


•  e:{p 


hit  + 


1  '  f  A  \  ^ 

— .  Bt  cos  $  sin  9.  >  , 


Y^(9,t)  =  •♦•  '—1  exp  {ift  + -—j  (-Bt  cos  6  cos  9)^/  ; 

^  \  z(ri  ri  23"^  1  ® 


t /p  /  T  \ -3/2  (Bt  cos  e  cos  9^) 

,(8,t)  =  - 5 - . 

\  2(rf  ^ 

111.  1  _  a  _  x2\. 

•  exp  \t;  t  (-Bt  cos^  cos  9  }  , 


(B12) 


(B13) 


(B14) 


i^ioltiplylng  Equations  (B12)  and  (BI3)..  and  (Bli)  and  (Blh)  together 
in  pairs  and  sicplifying  the  exponents,  the  follo^rtng  expressions 
axe  obtained: 

I  1  ^ A  \ 

x,(e,t)Y,(e,t)  =  -j-  t  ■  — jj  Bt  COE  e  sin  • 

^  ^  M  ze‘l  ''  «2l 


h  ^-r' I?  • 


COS  9  sin  9  +  B^t^cos^9 

^2  o 

0 


;B15) 
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r-2 


"t  +  '  *— jr 

2<r/ 


(Bt  cos  6  COB  9^) 


t  +  cos  e  SiB  (p  +  bVcO8^0 

26^  1(5*  o 


(B16) 


Kext,  in  the  eacponents  the  foUowiag  coi^lnatlon  is  seen: 


"M  .  /  .  1 


J 


-II 


.  (BI7) 


Thxia,  upon  substitutiisg  Equations  (BI5)  and  (BI6)  into  jSquation  (B3), 

and  waking  use  of  Equation  (BIT),  there  results  for  /e  iqpon 

2  2  \  o  / 

rearrangCTent  and  noting  that  sin  9^  +  cos  <p^  =  1 1 


/ 


dt 


exp 


1 


t  - 


I 


-II 


Bt  cos  &•> 


B^t^  I  1  '•  2  '  /  i  I 

•  cos  S  «aj  ^ cos  6  -  ^t|t  + 


26 

i 

where  again  the  paraaeter  P  is: 

fl  AB  ,  ^ 

p  "  • 

The  integration  with  respect  to  0  is  carried  out  in  terms 
of  the  following  basic  integral,  denoted  by  S(p,q); 


B 

(SIS) 
d0  , 


{BI9) 


( 


S(p,q)  =  J 


CJQ)  (p  COS  0  -  q  cos  0)  d0 


(B20) 


It  is  seen  from  Equation  (BI8)  that  the  desired  Integrals  will  be  obtained 
froa  Equation  (B20)  by  differentiation  with  respect  to  p  and  q.  The 

*  a 

evaIxMtion  of  Equation  (B2C ) ,  which  does  not  appear  obtainable  in  closed 
fora,  proceeds  by  expaadlDg  eacn  escponent  in  a  power  series: 
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exp(p  000^6) 


2  cos^^  e  , 

R=0 


(B21) 


e3tp(-q  cos  0) 


(B22) 


Puttixig  into  Equation  (B20)  and  interchanging  eke  order  of  opesrations^ 


there  results: 


«  00 


sC.PjQ)  =  21  2 


,  s»  n  m 


n=0  31=0 


b:  n; 


cos^^^0d0  .  f-23) 


Now,  breaking  up  the  integral  into  two  parts,  one  with  a 
range  0  -  x/2  and  the  other  x/2  -  x,  translating  0  in  the  latter  by 
-  x/2,  it  can  be  evaluated  by  Reference  19,  p.  196,  no,85^.1,  for 
m— )2n-Ha,  x  =  0,  yielding  upon  collecting  terms; 


d9  =  r^cos2"«  SiS^f 


0  d0 


fx/2  rs/2 

2n-Ha  o  ft  ,  Xx  ,« 

=  cos  0  CLW  +  cos  .(0  +  •?•;  d0 

J  ' 

>0  --  o 

^x/2  r  )t/2 

I  2n+ffift  ,ft  .  /  ^A2n-Ha  1  .  2n-Haft 

=  cos  0  d0  +  C-1}  I  sin  0  d0 


.Sn-Ha 


ml 

eae  =  Vf,  '  ,  ! 

i  Vn  +  1  +  ^) 


(BZk) 


But  it  is  seen  that 


1  +  (-1)“  =■  ( 


fo;  for  m  =  1,3>5>  . ^ 


2;  for  m  -  0,2,4,  . even  . 


(B25) 


Thus  setting  a  =  2k,  for  k  »  0,  1,  2,  ...  , 
Equation  {B24); 


there  results  instead  of 


I 
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2(n4k) 


e  dd  = 


r(n  +  k  +  i) 


[n  +  k  +  Ij 


(B26) 


Upon  substituting  Equations  (B24)  through  (B26)  into  Equation  (B23)^ 
it  is  seen  that  only  the  even  terns  in  the  sunmation  over  m  contribute. 
Thus,  setting  m  =  2k  here  as  well,  the  following  form  for  S(p,q)  is 


obtained: 


^n  2k  r(n  +  k  +  |) 


s(p,q)  =  --"-^^l:2^y:^r(nVkT-ij 


U  .  (B27) 


Next,  the  suosation  over  n  is  perforaed  first,  yielding: 


S(p,q) 


The  latter  series  is  recognized  in  terns  of  the  esqpaxislon  for  the 

( 

confluent  hypergeonetrlc  function  from  Reference  21,  II,  p.'  6, 
no*  (32),  for  a  »  k  +  I/2,  c  *  k  +  1,  z  =  p: 

"  ,  Rn  +  k  +  |)  Hk  +  |)  , 

^lTS-'k-i'i)  ^  •  rtrr4j •  (®'5) 


Putting  Equation  (B29)  into  Equation  (B26)  gives  the  expansion 


formula: 


»  ,2k  Rk  +  I) 


S(p„)  =  J,)  .  (B30) 


The  desired  integral  evaluations  in  Eq^iation  {BI8)  are  obtained 
from  Equation  (B30)  by  suitable  differentiation.  These  are: 


2  o 

exp  (p  cos  0  -  q  cos  9)  cos  9  d.0  =  -  ^  S(p,q)  ,  (B31) 


exp  (i>  cos  0  -  q  cos  0)  cos  0d0=  ^S(p,q)  .  (B32) 

<7^ 
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The  required  derivatives  of  Eqviation  (B30)  are  obtained  with  the 
aid  of  Reference  22,  6.4>  p.  254,  no.  (8),  for  a  =  k  +  1/2, 
c=k  +  l,  x  =  p,  *  stands  for  and  making  use  of  gassna  function 

reDjitions : 


^  S(p,q) 


=  n 


,  00  ^2k-l  Rk  +  0 

?2k-lji  p(k  +  1)  +  g;  k  +  1;  p) 


=  - 1  1=  -  2;  p)  .  (B33) 


«  2k+l 


Ri  + 1) 


=  It 


1/2 


®  2k  r(k  +  ^)  (k  +  i)  » 

ksTnOrriT  nrriT  - 1  “  -  2;  p) 


'  kri  -  fi  *  2;  p)  •  (B3M 


V,  «  „2k  Rk  +  I) 


Finally,  making  the  substitutions 

,  \-l 


P  = 


bV 


t  + 


;  q  =  Pt 


t  + 


2(5^1 


-1 


in  Equations  (B31)  througl-  (B34)  and  sxibstituting  appropriately  into 
Equation  (BI8)  and  arranging,  there  results  the  following  expansion 
formula  for 


00 

•  z 

k=0 


.2k 


T2k+1) I 


t 


1 

26" 


\-2k 


"  fi.  I  3 

pCk  +  2)  l^ip  ■"  ?' 


k  +  2; 


bV 


T 


t  + 


•  (2k  +  1)  B  + 


(B35) 
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Applying  the  gaana  fimction  relation  in  Reference  22,  1.2,  p.  Jio*  (^5) 
for  z  =  k  +  1,  it  is  found: 


(2k+l)J 


r[2(k  +  1)] 


22(k+l)-l 

"77^ 


Rk  +  1)  Rk  +  1  +  |) 


22k+l 

77^ 


ki 


Ric  +  |) 


(B36) 


Putting  into  Equation  (B35),  cancelling,  and  collecting  tenas,  noting 
that  P(k  +  2)  =  (k  +  l)i,  takes  the  fora: 


While  the  shove  form  has  heen  considered  for  approximate 
evaluations,  it  is  more  convenient  for  exact  numerical  integration  to 
make  the  change  of  varlahle: 


dx 


I 


26^(1  -  x) 


(B38) 


At  the  same  time  the  Kunmier  trazisformation  is  £q>pLied  to  the  confluent 
hypergeoaetric  function  from  Reference  22,  6.3,  p.  253,  no.  (7),  for 


a  =  k  +  3/2,  c  =  k  +  2,  X  =  (B^t^/k) 


for 


t  + 


and  ^  standing 
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A 


k  +  |;  k  +  2; 


=  exp 


bV 


t  + 


1  I 

is") 


-I 


It  * 

- j  ( 

1  20^1 

bV 


-i\ 


.  (B39) 


CoalDlnlng  the  exponents  containing  B  to  yield: 

i-l) 


-1 


and  then  applying  Equation  (B38)  to  Equations  (B39)  aad  (Bto)  vhen 
substituted  into  Equation  (B37),  the  new  expression  for  contains 

integrails  over  the  range  0  to  1: 


(Bl«) 


B 

2 


>0 


dx  €*^  Z 
k=0 

,2k 


■ 


(2k  +  1)  B  +  (1  .  x) 


•  kuA  1*  1^1  i*'i  (2- ^  t!:jJ 


where  the  parameter  a  is  defined  as  in  Appendix  A: 


^  1  f  a2  ^ 

a  =  — 7:  1 A  +  "5” 


2(5^ 

Equation  (b4i)  is  used  as  the  starting  point  for  the  evaluation  of 
^  quadrature  case  discussed  in  the  main  text  and  treated 
in  i^ipendlx  C. 


(Bl«) 


STL/!M-61-0000 -19008 


AFFE!iDI}C  C 

APPRODOMAIE  EVALUATIOSS  FOR  THE  OUTPUT  NOISE-TO-NOISE 
RATIO  FOR  THE  ^ADRATURE  CASE 


The  calciilAtiorx  of  froa  Equation  (27)  ty  nmaericeQ. 

integration  for  a  range  of  reference,  signal,  and  noise  paraneters 
becomes  quite  tioe  consuujjog,  so  it  is  desirable  to  seek  accurate 
approxiaate  evaluations  vhere  possible.  The  derivations  presented 
here  cover  the  cases  of  saall  and  large  reference  ac^litudes. 


Snail  Reference  Aaplltude 


ExaaixiAtion  of  Equation  (27)  for  y  would  indicate  at 
first  look  that  for  B  <  <  2^^^6,  it  is  aerely  enough  t<o  take  the 
first  few  terms  in  the  series  e:q)ansion  for  the  confluent  hypergeonetric 
function,  found  in  Reference  22,  6»1,  p.  2^,  Ho.  (l)  for  a  =  l/2, 
c  *  2,  X  replaced  by  -  (bV8(5^)(xVi*3c): 


1^ 


1.  2- 

-  1  Fb^  x^ 

^  1 

B^  X^ 

2'  '  '  88^  TI3J 

^15 

[88^ 

■  1^5  [ 


sL 

88^ 


(Cl) 


But  an  Integration  between  the  Halts  0  and  1  mxist  subsequently  be 
carried  out,  and  the  above  series  obviously  dl. verges  as  x — )1. 

However,  closer  exaal nation  of  Equation  (Cl)  reveals  that 
the  square  of  the  second  tern  gives  the  third  term.  With  the  above 
series  alternating  in  sign,  and  denoting  by  y  the  quantity: 


y 


} 


(C2) 


vhere  t  is  the  normalized  ratio: 


0X1*/  xn*»ox-vwv  -x^^vX>u 


Page  2C 


r 


(C3) 


then  it  Is  seen  that  the  first  three  terms  of  Eqxiation  (Cl)  are 
identiceJ.  with  coirespondlng  ones  of  the  alternating  series  found 
in  Reference  19>  p.  3^  Jio*  9*0^^  y  =  x; 


1 


y  -»■  y 


1 

1  +  y 


ick) 


Further,  the  fovirth  tem  of  Equation  (Ck) ,  which  in  terms  of  x 
from  Equations  (C3)  and  (c4)  becomes 


> 


differs  from  the  corresponding  term  Of  Equation  (Cl)  by  about 
13^,  idiich  for  smeQJ.  y  would  be  negligible  over  a  considerable 
range  of  y«  Ignoring  for  the  moment  'tiiat  the  series  in  Equation 
(d4)  does  not  converge  for  y  ^  1,  it  would  appear  that  a  possible 
approximation  to  for  .small  r  is: 


1-x 


I 


(C5) 


It  still  remains  to  see  how  the  approximation  behaves  in 
the  range  1  =  y  =  ••  corresponding  to  values  of  x  approaching 
unity.  At  X  =  1,  the  denominator  of  the  right  side  of  Equation  (C3) 
becomes  infinite,  so  the  over-all  expression  is  zero.  From 
Reference  22,  6.13»lt  P»  278,  no. (3),  it  is  seen  that  as  the 
argument  of  ^F^  approaches  -  «•  as  x — >1,  the  8U5ya5)totic  behavior 
of  the  confluent  hypergeonetric  function  for  the  parameters  l/2. 


2  is: 
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1 

1?  i*  ?• 

r  1  z* 


«2\  9 

-  el  h  -I 


— t  0  as  X  —  ■'  >  1  • 


Hms,  the  original  runction  aad  the  proposed  approxlaatloa  agree 
In  this  jjQ>ortai)t  respect. 

To  concli^lvely  check  the  asciuracy  of  2q[uatlon  (C^)^  It 
Is  plotted  against  the  actual  function  for  two  ralues  of  0.1  and  1, 

on  Plgure  Cl,  over  the  range  0  «  x  «  1.  The  case  r  “  1  i® 
considered  i^pper  Unit  for  the  approximation,  for  vhlch 
Equation  (C^)  reduces  to  the  special  form: 


I'Jl'  ^ 

^  ^  *  ^-*1  (Z-xf 

The  exact  valties  of  were  obtained  with  the  aid  of  tables  In 
Beference  23#  PP«  698-TI3#  since  it  can  alternatively  be  eacpreesed 
in  terns  of  esq^nentlals  and  modified  Bessel  functions  by: 


,-t/3 


yi 


+  I. 


which  can  be  derived  from  relations  given  in  Reference  22,  Chapter  6. 
The  very  close  fit  for  the  x  -  0,1  case  plus  the  good  agreement  for 
the  limiting  case  gives  confidence  that  the  approximation  will  be  quite 
accurate  aver  the  range  .1^1  1#  aad  even  better  for  smaller  r* 

Equation  (C5)  is  applied  to  the  evaluation  of 
Equation  (2?)  by  rewriting  the  right-hand  side  in  xiormal  ratio  form: 


rx^  -  %  +  h 


fl-x) 


The  roots  for  the  denominator  quadratic  in  x  are: 
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vhlch  are  real  and  distinct  for  r  ^  with  equality  considered  as  a 
special  case  hy  Equation  (C6).  lEhe  above  restriction  for  real  roots 
thus  coincides  with  the  previously  considered  iqjper  Halt  for  the 
aK>roxlaation»  How,  the  right  side  of  Equation  (C8)  can  be  written 
‘>lth  the  factored  denoainatcr,  tising  Equation  (C9): 


1  .  4^1-x) 

^  ..2  ‘  r(x-h.)(x-h  ) 


(CIO) 


This  expression  can  be  expanded  further  by  partial  fractions;  which 
after  some  coaiputation  and  canbinlog  with  Eqiiatlon  (C5)  yiej.ds  the 
desxTed  fora  of  approximation: 


ri 


(i.  2-  - 

-w  r 

[2  2-1 

(2^  (l-x) 

4(i»r)^^ 

L(x-u^)  '  (x-u^)J 

(cu) 


Vpon  stibttitutlon  of  Squotiou  (CXl)  into  Equation  (2T)  with 

/  2\ 

aquation  (C3)  for  the  pcjtRBster  y,  there  results  for  ' 


arrangement : 


_ A- 

8(1?)^ 


!  zf 

i 


1 


2 


e  ")  upon 
®  / 


6  dx 


,'0 


(012) 


Since  for  0  r  n  1  the  inequality  >  i  bolds,  the  denoalnstois  of 
the  Integrands  In  Equation  (C12)  do  xiot  vanish  over  the  range  of 
integration,  and  thus  the  Integrals  exist.  They  can  be  evaluatf^d  in 
terms  of  the  modified  expcnenllal  integral  £i(x),  defined  in  Beference  12, 
9»7f  P«  1^3,  iio.(3)^  as  a  Cauchy  principal  value: 


KICa) 


.-t 


dt 


.€ 


, ,  r  I 

i.e.,  lltt  (  /  t  I  / 

j.  Li 


-X 


as  <5 — >0  * 


'>-r 


(013) 
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By  making  the  translational  changes  of  variabj.e 
t  =  a(x  - 

in  the  appropriate  integrals  of  Equation  (C12'j  breaking  them 
up  into  two  coaqwnents  with  infinite  upper  limits  an.!  applying 
Equation  (CI3),  there  results  upon  arrangement  the  evaluation; 


--ax, 

■  ■  -  <-.) 


Jf  ‘-f-f 


=  -  exp  (-04^)  Ei(a4^)  -  Ei  Mu^-l) 


(cik) 


Substitution  of  Eq\»  “Ion  (Cl4)  into  Equation  (C12)  yields  the 
expression  for  )  fo-  0  ^  r  ^1* 

Foa  the  limiting  case  r  ”  the  approximation  to 
given  by  Eqvustion  (C6),  which  takes  the  expanded  form; 


1  x^  1 

IP  i  •  2  _ _  =’  4 

f  1  2  ^  ^  ^  1-x) 


U-2) 


1  1 

’ 


(CI5) 


Putting  into  Equation  (27)  gives  for  ^^^^in  this  special  case: 


/e  =  2B^ 

\  o  / 1 


;  -ax 


f 

.2)^'  ■  I 


-a.x 

€  ax 


(x-2) 


(Ci6) 


where  integral®  agsdn  exist.  Integrating  by 

;  ^i'^s  in  the  first  term,  with 


u  =  € 


dv  = 


(x-2)‘ 


yields  upon  coBibining  terms  the  form: 
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-  (l  + 


“‘1 


-CSX 
€  dx 


Ti=zT 


(C17) 


!Ihe  latter  integral  is  evaluated  again  In  taros  of  aodified 
e:Q>onential  integrals  as  before  by  the  cbaage  of  variable  t  =  (X^Cx-s) 
to  give; 


-V  -**i 

€  dx  ^ 

■» 

f  e'*  ^  f 

1  ^ 

1 

—  »*■  ~f  —  « 

Jo 

LJ 

-ac^  J 

-aa. 


-  € 


Bi(2a^)  -  Ei(a^) 


(Cl8) 


Substituting  Equation  (ClB)  into  Equation  (C17)  gives  for  ' 


+  € 


Z!  ,  (C19) 


vhere  frwa  the  definition  of  a  by  Equation  (21)  and  r  by  Eq^uation 
(C3): 

=  1  +  ^  .  {C20) 

To  coa^lete  the  calculrtion  of  the  noise -to-noise  rati’>, 
it  is  necessary  to  evaluate  the  ej^ression  for  the  aean  square 
value  for  aero  signal  iwp^^bude^  A  =  0.  Ihider  this  condition^  he 
foUoving  paraaeter  reductions  occur,  denoted  by  zero  stibscripf : 


%  =  r  !  »  1  .  (021) 

Putting  appropriately  into  Equations  (C12)  and  (Cl4),  the  fonsula 
for  at  A  =  0  is  obtained.  Finally,  the  substitutions  into 

Squjitionci  (26)  are  aade  vltt  s  o  to  yield  the  £i)proxiBate 

cjqpresslon  for  5  /ff_  • 


d 
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Saall  Befereace  Amplitude  Llait 

Here,  the  appropriate  foraaila  for  found  from 

the  limiting  fora  found  in  Appendix  D,  Equation  (D30),  for  the 
general  ce^e,  using  Equation  (30)  and  noting  from  Equation  (2l) 
that  as  B— »0  : 

a  — )  ^  =  a  .  (CE2) 

2(r 

Making  the  appropriate  substitution  in  the  above-mentioned  expression 
and  setting  (p  =0,  there  results  the  expression  appearing  in 
Equation  (33)  of  the  main  text  for  since  as  A — >0,  the 

iade terminate  evalmtion  >  elds  the  factor  B^/2 : 

.  (023) 


Leuge  Beference  Aag?itude 

Ihe  use  of  Equation  (27)  as  a  staurfcing  point  for  attempting 
an  approximate  evaluation  of  large  B  led  to  difficulty 

in  assessing  the  validity  of  the  assm^rtlons  made.  Further ^  a 
reduction  of  the  corresponding  limiting  formula  in  Appendix  D 
for  <p^  =  0  gave  an  oversin5)lified  result  of  no  value.  Thus^  to 
obtain  \iseful  e:q)ression8,  an  asynqrtotic  expansion  of  the  original 
I^iase  detector  output  voltage  was  taken  as  the  basis  ^r  the  following 
derivation. 


Noting  from  Equations  (8)  and  (9)  that  the  voltages  e^  mid 
eg  from  the  circuit  halves  differ  only  by  a  sign  in  one  term,  they 
can  be  rewritlien  in  combined  fora  as  e.  after  factoring  out  B/2: 


e 


1 


o 


B 

2 


l  +  sln 


(C.24) 


\ 
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It  ia  nov  assumed  that  the  reference  aa5)litudB  IJ  is  sufficiently 
larger  than  the  input  sigxial  plus  noise  asqclitude  p  for  all 
values  of  9  to  aahe  the  infinite  series  e3q)ansion  of  (l  + 
valid^  vhere: 


2  =  jf  ±  sin  <p 


(C25) 


Fro*  Beference  19>  p*2,  no.  5*3>  ^or  x  =  z,  the  condition  is 


z^  =  1 


and  the  expansion  is: 


(1  +  z)^^  =  1  +  i 


^  1.1  2  l-l-; 

^  “  275  ^ 


X'  - 


l-l.; 


k 

z  + 


(C26) 


Reference  Zh,  p.ll7,  no.  6.33,  for  z  =  x  gives  the  general  tera  in 
this  e3q>an8ion  for  k  «  2: 


(-1) 


k  * 

k.’  2^ 


But  since 


1.3. 5....  (2k-3)  =  -r-  , 

Equation  (C26)  can  be  revritten  in  suanation  fona  as: 

.  l.iz  -  Z  _ 

k«2  .2^"\:{k-2)! 


(C27) 


Raising  Equation  (C25)  to  the  k-th  power  and  applying  the 
binominal  expansion  of  Reference  19,  p.  1,  nos.  1  and  for 
a  =  p/b,  x  =  sin  9  and  n  =  L: 


2  ;  for  m  odd. 
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Setting  a  =  2n  +  1,  so  thet  for 

«  =  ;  a  =  [^]  ,, 

where  the  i^tter  notation  laeans  that  if  k  is  odd^  it  is  the 
upper  limit  of  - the  sua^  hut  if  k  is  even,  then  the  sum  terminates 
at  the  next  lower  Integer  index,  (k-2)/2.  Thus,  in  place  of 
E4uation  (C30),  there  results  with  the  index  n: 


2p  sin  (p 


2  It  (-l)^(21c-3)!  I 


,2n-fl 


2(k-n)-l  ._2n-M^  _ _ _ 

riy;(k-2n-i)!  0  ^ 


Writing  out  the  first  few  texms  of  the  expansion  yields: 


2p  sin  9  -  3  ^  l!n  ^ 

V  '  I 

*3^  (iiIt ^  ^ 


.  iL  „7 


B  5.3 


T2;  1T3TP  9  +  37X1  P"  «?] 


(C32) 


^  T  *  B  5  B^  T  ^  B  ^  S  S 

1>  +  ^  pWip 


•  If  (nsT  ^  ^ 

^  ^  (i75T  ■■■  ■'■  i4t  P^Pi^’^l'l 
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Multiplying  out  the  numerical  factors  and  collecting  terms  in 
inverse  powers  of  B  up  t..  the  following  expression  for 

e^  is  obteiined  after  arrangement: 

e  =  2p  sin  <p  -  (p^  sin  9  -  p^  sin^  Cp) 

o 

+  \  (Tp^  sin^  9  -  lOp^  sin^  9  +  Sp^  sin  9)  (C33) 

B 

+  ^  (33p'^sin^  9  -  63p'^sin^  9  +  35p^sin^  9  -  5p^sin  9)  +  ••• 
B 


In  order  to  con^jute  the  output  noise -to -noise  ratio  for  the 

quadrature  case,  only  the  ensenOble  average  of  the  square  of  the 

output,  (e  ^),  is  required.  So,  i;q>on  squaring  Equation  (C33) 

and  collecting  terms  up  to  1/B  only  for  the  asyaqrbotic  evaluation 

2 

considered  here,  the  approximate  expression  for  becoases: 


^  4 


^  sin^  9  -  (p^  sia^  9  -  P^  sin^  9) 


^16  /^6  _6  6  .2^\ 

+  -r  {Zp  sin  9  -  3p  sin  9  +  p  sin  9) 

B 


} 


which  furthe.  t.’'i''onoMetric  sinqjlification  reduces  to  a  form  more 
convenient  i«  averaging.  Noting  also  from  Equation  (lO)  tiiat  for 
the  qviAdrature  -'ase,  9  =  9^,  it  is  found: 


Cl  ^  \. 

e  =  4 

o 


2.2^  ^/2.2^w2  2  ^  \ 

p  sin  9^ - ^  (p  sin  9^)Cp  cos  9^) 

B 


+  ^  p^  sin^  9^(2  sin^  ~  ^ 
B 


(C34) 


The  ensemble  averaging  of  Equation  (C34)  is  sija5)lified 
considerably  by  the  fact  that  some  of  the  factors  can  be  expressed 
directly  in  terms  of  the  in-phase  and  quadrature  noise  conqxonents 
and  the  signal  amplitude  by  referring  to  Equations  (2)  and  (4). 
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From  them  the  folloving  relations  are  iouad: 


p*^sin^(p^  =  ;  (p"^8in^?^)(p‘^coB^q)^)  =  y^(A+x)^  ; 

p^8ln^9^  ® 

Making  the  appropriate  aubstittxtions  in  Equation  (C3^)  and  taking 
the  average  on  both  sides,  there  results  for 


S'  4 


(C35) 


{/)  -  ^  (C36) 

^  L  -  3  (p^B±n\^  +  ^^sin^(p^\|  . 


B 

The  terms  expressed  in  rectangular  conponents  are  treated 
first*  Since  x  and  y  are  Independent  Gaussian  random  variables 
with  mean  value  zero  and  standard  deviation  6,  ther.  v'^th  the 
aid  of  Jteference  25>  p*  83#  no.  10. 5>  the  following  u  jrages  are 
fotuxd: 


(C37) 


<x)  =  0  j  (x^)  =  0;  (x^)  =  0  ;  (y)  =  0  ; 

(x^)  =  ;  (y^)  =  ;  (y^  = 

Further,  since  A  is  a  constant,  it  is  clear  that: 

(a^)  «  A^  .  (C38) 

CoflQnzting  the  average  indicated  in  Equation  (C36)  with  the  aid  of 
Equations  (C37)  and  (038)  yields; 

=  <yV)  .  2  .  (yV> 

=  (y^)  *  2A  (y^)(x)  +  (y^/(x^)  or 

aV  +  (5^  =  e^(A^  +  6^)  .  (C39) 


(y^{A+x) 


2\ 

/ 
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The  remaining  two  averages  to  he  found  in  Equation  (C36) 
must  be  determined  in  a  probability  sense  from  the 

JoiJit  density  function  of  amplitude  and  phase  of  the  signal  plus 
noise  given  by  Equation  (l3)*  resulting  double  integrals 

are  therefore: 


sin 


’i/ 


dp*p  I  sin  W^(p,<p^)  dfp^  , 


'O 


(p^Btri^(p^  =  I  dp‘p^  f  sin%^  W2{p,(p^)  djp^ 

"-O  -yQ 

Putting  in  Equation  (13)  and  arranging  yields: 

r" 

'6  .  k_\  exp  (-A^/2cr)  .  7 _  /  p^ 

p  c.xji  9.  )  a  ■■  *  ^  ^ ^  <3p  •  P  exp  (- 

1/  2x6^  ' 


f’  U 

Ap  CC3  9. 

sin  9^  exp 

) 

-'O  , 


/6  .  2^\  exp  (-/w2(5^)  [  ,  7 

(p  s.in  9^  /  =  — -i-i — - i-  !  dp  •  p  exp 

\  2ncr  .0 


80 


-  n 


r 

2 

•  2  ,  sin  9^  exp 
^  o 


lAp  cos  9. 


6" 


1  26^1 

d9, 


26^ 


d9. 


(ClK)) 


(Ci+l) 


Applying  Reference  20,  p.  202,  no.  175^  for  V  =  1,2,  z  =  Ap/(J^  , 
there  res^ilts  upon  cancelling  terms: 

CO 

/6  ,  i^_\  3<5^exp(-A^/2(5^)  !  5  I 

>(p  sin  9^j  =  ^ ^  [  p-^  ejqp  j- 

^  Jo 


!-  ^ 

I 

fie 

2 

dp  ,  (042) 


(A 


6  ^  2_.  \  expi~fi?/2C}^)  f  6  ^ 

p  sin  9j^^  =  — ^  '  p  exp 

J  o 


2  1 

.  E_ 

i  I 

26^1 

'  Ms" 

dp 


,  (c43) 


>rtiere  and  are  the  modified  Bessel  functions  of  first  kind  of 
first  and  second  ordtjr,  I'espectively. 
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*  change  of  variable  1;  =  p  presex'vljig  "the  llal'ts 

results  In: 


/S  ,  k  \ 

(p  sin  9^) 


-'o 


2A‘ 


*  tV2| 


exp  - 


(pW9^')  »  25EI^<£1  f  t5/2  1^ 


exp  - 


V 

2(ri 


dt 


at 


Integrals  of  the  above  type  can  be  evaluated  generally  in  terms 
of  confluent  hypergeoaetrlc  functions  by  Reference  9,  table  4.l6, 

P»  197;  BO.  20,  and  then  reducing  the  results  for  special  parameter 
cases.  But  a  more  direct  determination  is  seen  by  suitid>le 
differentiation  of  220  •  18  of  the  same  table  4>r  V  =  1,  2 »  This 
gives  the  foUovlng  fomulas: 


dt 


t  I-(aa^/V/^)€"P^dt 


«  A 


(fj 

=  ^  (3p-kx)  exp  j~j 


dp  \p 


P// 


0^2 


(C44) 


(c45) 


{Ck6) 


(6p^+  6ttp  +  c?)  exp  (2j  , 


Setting  the  parameter  values 


1 

2<5^ 


o 


STL/TM-61-0000-19008 
P«"S  2j6c 


In  Equations  (C46)  and  (C47)>  sioipllfylDg,  and  then  suhstituting 
appropriately  into  Equations  (C44)  and  (Ck^),  there  results  upon 
cancelling  terms  the  ensemble  averages: 


(p^sinS^ 

(p^sin^cp^) 


j 


6  +  6 

0 

+ 

>»> 

to 

■ 

26 

26*^'  J 

(ckS) 

(c49) 


Finally,  applying  Equations  (C37),  {C39),  (Chd),  and  (C49) 
to  Equation  (C36)  and  arranging,  it  is  found  for  the  s»an  square 
value: 


Ka2d.ng  use  of  Eqizations  (28)  and  (30)  wheu  the  dioensionless 

ratio  V  and  the  input  nignal-to-noise  power  ratio  "a'’  are  defined 

,  2\ 

respectively,  the  fonnula  for  }  for  large  reference  amplitudes 
takes  the  simplified  form; 

('o'>  =  “s' 

From  this,  the  formula  at  A  =  0  is  obtained  by  setting 

a  =  0,  '.diich  when  substituted  together  with  the  original  expression 
into  Equation  (26)  along  with  (&\  =0  yields  the  approximate 

expression  for  5  /6  appearing  as  Equation  (3^)  of  the  main  text. 

9 


ii  - 


l+2a  ^  1 
2r 


2r^  ^ 
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ATJtJPimTV  n 

*  haaia#  j  #».  a/ 

LDOTIHG  FORMULAS  FOR  THE  OUTPUT  SIONAL-TO-NOISS 
RATIO  FOR  TEE  GENERAL  CASE 


The  approxlaate  evaluations  of  Equation  (B37)  considered 
for  resulted  from  a  ntaaber  of  coaqmtational  steps  vhich 

tended  to  obscure  the  nature  of  the  approxiaations  made  and 
thus  cast  doubt  on  the  range  of  their  validity.  Therefore, 
it  appeared  desirable  to  find  some  alternative  derivation  to 
obtain  useful  foraxzlas.  A  very  direct  method  presented  itself 
in  vhich  the  starting  point.  Instead  of  involving  t^ie  integrated 
results,  is  the  original  phase  detector  output  voltage  expression. 
Taking  (8)  and  (9)  as  the  voltages  from  the  circuit 

halves,  squaring,  <*ubtraetiag,  factoring,  and  using  Eqiiatlcxas  (lO) 
aad'(ll),  there  results  the  following  form  for  e^: 


1*1  +  e^j 


sin  (q>^  + 


From  this  relation,  two  limiting  cases  in  terms  of  the  reference 
SBiplitude  are  derived  for  the  output  signal-to-noise  ratio 

Large  Reference  Amplitude  Limit 

If  the  reference  asplltude  B  is  38uch  larger  than  the 
aaqpUtuie  p  of  slgxial  plus  noise,  then  with  the  aid  of  Figure  2 
it  is  seen  ’that: 


*1  *  *2  "  ® 


!Qxus,  for  the  limiting  ease.  Equation  (D1)  reduces  tc: 


2p  sin  (9^  +  9^) 


independent  of  B.  Escpandlng  and  grouping  terms  yields: 
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2[{p  sia  cp^)co8  <p^  +  (p  cos  q>^)  sin  (p^J 


But  the  quantities  in  parentheses  can  be  expressed  directly  in 
terns  of  the  in-phase  and  quadrature  noise  components  and  the 
signal  amplitude  by  referring  to  Eqxiationa  (2)  and  (4),  which 
results  in  the  following  simplified  fora  for  e^: 

e^  =  2  [y  cos  (p^  +  (A+x)  sin  <p^j  .  (D5) 

{ 

Upon  taking  thC  ensemble  average  of  Equation  (D^)  and  its 
square,  noting  that  factors  depending  only  on  <p^  are  Independent  of 
this  operation,  it  is  found: 

(*o)  “  ^  [<^)  '^o  (^)  %]  * 

cos^(p^+  2^<^)  +  sia  (p^cos  (p^ 

+  ^  sln^  9^]  .  (D7) 

Fzx>a  the  facts  that  x  and  y  are  independent  Gaussian  random  variables 
with  mean  value  zero  trad  etandard  deviation  6,  and  that  A  is  a 
constant,  the  following  relations  hold: 


(*>  =  0  ;  <;^y)  -  0  ;  =  0  ; 

■■  0  i  »  0  ; 

(a>  =  A  j  <A^  =  a''  ;  =  8^  ;  (/}  ^  (^ 

Substituting  into  Equations  (d6)  and  (D7)  yields: 

<«o)  =  2A  sin  9^  ; 

i  k{(^  4  A^  sin^^) 


(D8) 

(D9) 

(DK») 
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Applying  the^-j  'osultfl  to  Equation  (36),  siiiplifylng  and 
■akiag  use  of  Equation  (^)  ^ere  the  input  signal-to-noise  power 
ratio  "a”  is  defined,  -^^ST  is  found  to  have  the  limiting  formula 
in  the  large  reference  aji^litude  case; 


8in(p^  .  (Pll) 

This  is  found  as  Equation  (37)  of  the  main  text. 

Small.  Reference  Amplitude  Limit 

If  the  reference  asQ>litude  B  is  much  smaller  than  the  aaq>litude 
P  of  signal  plus  noise,  then  with  the  aid  of  Pigm-e  2  it  is  seen  that; 


*1  *2  “ 

Thus  for  this  limiting  case.  Equation  (Dl)  reduces  to; 


(DI2) 


«o  *  S  •ia  +  %)  , 

independent  of  p.  Upon  ejQWHadlag  Equation  (DI3); 


(DI3) 


Cq  »  B(sin  9^  cos  9^  +  cos  9^  sin  9^)  ,  (D14) 

it  is  seen  that  the  factors  involved  in  the  statistical  a»«raging 
cannot  he  reduced  into  sii^jle  rectangular  forms  as  before.  So, 

carrying  out  the  averaging  on  Equation  (PJ4)  and  its  square, 
it  is  found: 


cos  +  ^OS  tlB  9oj  , 

Tq  *  ^  fi*®*  <P^ 

•ij?  9^1  . 


(DI5) 

Sin  9  cos  9 

o 

(PI6) 
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Asaired  Averages  are  found  from  fommlaa  analogoiis  to 
Equitiona  (l4)  and  (15);  vhcre  for  a  general,  function  of  «p^, 

2X 

(D17) 


^  2k 

(tm)  -  I  a,o  I  f((p^)  d/p^  , 


where  V^(Pf<p^/f  the  Joint  probability  distribution  of  aiaplitude 
and  phase  of  the  signal  plus  noise  is  given  by  Equation  (13)> 

Its  application  to  Equation  (Dl?)  for  the  trigonometric  functions 
indicated  in  Equations  (DI5)  and  (D16)  follows  in  order  of  their 
appearance . 


00 


dp  .  p  exp 


2*6^ 


^in  <p^  *=  0 


r2K 

lAp  COB  9. 

1  sin  9^  exp 
'0 

[  ^ 

1  ' 

-  jlA. 
2^1 


djp^  or 


(DIB) 


since  the  integrand  in  is  an  odd  function  over  the  0  -  2jc  interval. 


.  exg(-^/2a^)  j  ^  . 


P  exp 


r  2j( 


cos  exp 


Ap  cos  <p. 


2 

-  £— 
2(5^ 


dcp. 


Applying  Reference  20,  p.  202,  no.  I80,  for  n  =  1,  z  -  Ap/6^ 


<00.  9,)  =  [  p  J 


Jo 


exp 


2 

.  £_ 

2(5^  ^ 


clp  ,  (DI9) 


where  I,  is  the  modified  Bessel  function  c.f  first  kind  and  order. 

2 

Making  the  change  of  variable  x  =  p  in  Equation  (DI9}  results 
in  the  form  with  the  sane  limits : 
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/cob  9.\ 
\  -1/ 


exp(-A^/2l5^)  f 

-'O 


{D20) 


Thl#  integral  io  evaluated  in  Reference  l8,  p.  63^  no.  3«^31.»  for 
p  »  1/2(5^,  A  *  ^ji/ ,  V  =  1,  a  »  3/2,  which  upon  substitution  into 
Equation  (D20)  and  sia5>llficetion  yields  for  <co8  cp^  ; 


\1tA  _ 


(D21) 


where  is  again  the  go/ fluent  hypergeoa»tric  function.  For 
the  particular  parajaeter  "^ound  above;  this  function  can 

also  be  expressed  in  tcxffis  of  :,xp-inentials  and  aodlfled  Bessel 
functions  fx«a  relations  given  in  Reference  22;  Chapter  6;  as 
was  the  case  in  Appendix  C,  Equation  (C7).  The  resxiit  is: 


F 

ri 


2*  t 

•3*  6,  ^ 


gt/2  iti  .  ^  iti 


(D22) 


Substituting  into  Equation  (D2l)  with  t  =  A^/2(5^  th”.\s  g'.ves: 


’1)  ■  ^  $ 


+  I, 


which  is  a  sore  convenient  fora  for  csCculatlon. 

00 


<■ 


2«(5‘ 

r 


■Jo 


dp  .  p  exp  I  -  ^  I  • 
25  j 


|Ap  cos  f. ^ 

2  /  sin  vi)^  exp  ; - 5 - ,  djp^ 

^o 


(D23) 


Applying  Reference  20,  p.  202,  no.  175#  for  V  *:  1,  2  =  Ap/cS^,  there 
results  upon  canceling  terms: 


( 


eaqp 


''O 


28^ 


^  • 

fyt 


inzk) 
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The  change  ot  variable  t  =  p  preserving  the  Umlta  yields! 


^  1 


exp  -  i  ^  .  (D25) 


^/?73 


This  Integral  is  evaluated  in  Reference  9>  table  U.l6,  p.  197^  no.  l6 
for  p  =  l/2CS^,  a  «  vhich  when  p\rt  into  Equation  (D25)  gives; 


(slnSj) 


exp^«A 

2A 


zi’  f  |A^  I 


=  4  [i  - 


ea<p  - 


a  form  convenient  for  later  parameter  change. 


(D26) 


<^in  qp^cos  = 


exp(-A  /2< 


dp  •  p  exp  j-  j  . 


^  wW  /  t,  . .  I 

Ap  cos  cp.  I 

•  sin  cos  (p.  esq?  - = — =1  dJp. 

11  ^  !  1 


^in  ^iCos  9^  =  0  , 

since  the  integrand  in  9^  is  an  odd  function  over  the  0  -  2x  interval. 

Finally j  the  remaining  ensemble  avlerage,  ^os^qp.),  can  be 

/  2  \  ' 

found  directly  from  the  result  for  <^sin  9^y  by  trigonometric  identity 
and  the  application  of  Equation  (D26): 

^  -  sin^9^^  =  1  -  ^in^9^)  or 

<^oos\)  =  1  .  ^  |l  .  exp  ^1]  . 


(D27) 


(D28) 
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Svbstitutiag  EquatioM  (Dio),  (ii23)>  and  (d26)  to  (D28)  appropriately 
Into  E«uations  (DI5)  ssd  (Sl6),  arrauglDg  and  Biji5)lixying,  there 
resultf  the  desired  enseaiblo  averages: 


+  I, 


sin  9^  , 


cos  29+  sin  <p^ 


(D29) 


(D30) 


Applying  these  results  to  Equation  (36)  and  making  iise  of 
Equation  (30)  where  the  input  signal-to-noise  power  ratio  ”a"  is 
defined^  the  Uaiting  formula  for  -y/aT  in  the  small  referexice  asqplitude 
ease  is  obtained: 


(D31) 


This  is  found  as  Equation  (38)  of  the  main  text. 


